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Abstract 

We study the consequences of supersymmetry breaking in the computation of the number of solutions 
of the Thouless- Anderson- Palmer (TAP) equations. We show that Kurchan argument that proves the 
vanishing of the prefactor of the Bray and Moore saddle point for the total number of solutions can be 
extended to solutions at any given free energy. We also provide a new simple argument for the vanishing 
of the prefactor and use it to prove that the isolated eigenvalue recently considered by Aspelmeier, Bray 
and Moore is exactly zero in the BM theory because of supersymmetry breaking. The behavior of the 
eigenvector of the isolated eigenvalue at the lower band edge is also considered. 

1 Introduction 

The Complexity, defined as the logarithm of the number of metastable states divided by the size of the 
system N, is a key concept in disordered system theory. Recently, the classical Bray and Moore (BM) 
computation of the complexity in the Sherringhton-Kirkpatrick (SK) model [7J has been reconsidered 
El I1U1 ITU] starting from the observation that it violates a supersymmetry (SUSY) as noted in J3| . In this 
work we discuss some consequences of SUSY breaking. 

Stable and metastable states of the system can be associated to solutions of the Thouless-Anderson- 
Palmcr (TAP) equations [2] but in general not any TAP solution can be actually identified with a thermo- 
dynamical state, see e.g. discussion in this point however is not important to the present discussion 

therefore for the sake if simplicity we will identify the complexity as the logarithm of the total number of 
TAP solutions. Much more important is the fact that for technical reasons the actual object we consider is 
not the sum of all solutions of the TAP equations but rather a sum where each solution is weighted with the 
sign of determinant of the Hessian of the TAP free energy. This modification is very important and it is the 
very origin of the SUSY of the problem. The first consequence of this modification is that the weighted sum 
over all solutions must be trivially equal to one, indeed the Morse Theorem states that the total number 
of solutions of the TAP equations with positive determinant minus the number of solutions with negative 
determinant is equal to one. Therefore the problem splits in two: from one hand we compute the sum 
without the modulus but on the other hand we want to extract from the very same computation the actual 
complexity, i.e. the sum with the modulus. In other words we want to compute the total complexity but 
we want also to recover the trivial result that the sum without the modulus of all solutions be equal to 1. 

In the following we briefly outline the results of the paper. In section [5] we will express the number of 
solutions of the TAP equations as an integral over the exponential of an action depending on bosonic and 
fcrmionic (Grassmanian) variables. Due to the mean-field nature of the model this integral can be evaluated 
for large N by the saddle point method. Then the prefactor of the leading exponential contribution can be 
obtained as a power series in 1/N. The problem is to solve the saddle point equations obtained in order to 
extremize the action. A solution of the saddle point equations was obtained more than twenty years ago 
by Bray and Moore (BM) [7J and another solution was proposed in The two solutions lead to a very 
different complexity, see e.g. |!|], and one must decide which is the correct one. The action is invariant 
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under the so called Becchi-Rouet-Stora-Tyutin (BRST) transformation |11U12| . This transformation mixes 
bosonic and fermionic variables i. e. it is a SUSY transformation. An important difference between the BM 
solution and the second solution is that the BM solution breaks the SUSY of the action, as firstly pointed 
out by Kurchan in while the second solution preserves the SUSY 0. 

Recently it has been shown ^Oj that the SUSY solution is inconsistent (except at the lower-band edge, 
see in full replica-symmetry- breaking (FRSB) models like the SK model (at variance with 1RSB 

models), thus the BM solution remains the only viable candidate to describe a finite complexity in these 
models and the knowledge of the consequences of SUSY breaking becomes very important. 

As a first consequence of BRST invariance, one finds that the the prefactor of the exponential contri- 
bution to the complexity vanishes when evaluated at a non-SUSY saddle point, like the BM one. This 
result was obtained by Kurchan in |13j for the total number of solutions, here we further extend it to 
solutions at any given value of the free energy. The result is obtained in two steps. In the next two sections 
we start recalling that the expression of £(/) (the logarithm of the number of TAP solutions of a given 
free energy /) is invariant under a microscopic supersymmetric transformation that is a generalization of 
the standard BRST transformation [^j . After performing the disorder averages the integral over the micro- 
scopic action is expressed as an integral over a macroscopic action depending on eight macroscopic variables 
(four fermionic and four bosonic), that are introduced through standard manipulation such as even and 
odd Hubbard-Stratonovich transformations. Following |13| we note that this action too is invariant under 
a macroscopic SUSY transformation (that depends on the parameter u conjugated to the free energy /) 
between its eight variables. The second step consists in applying a general result that states that if the 
action is invariant under a SUSY transformation then the prefactor of the exponential contribution is zero 
at all orders in an expansion in powers of 1/N. This result is completely general and was obtained by 
Kurchan 13 who applied it to the total number of solutions. Since we have shown that the macroscopic 
action at any give value of the free energy is also BRST invariant we can apply it in this case too. For 
completeness Kurchan argument is reported in section 03 

In section0]we provide a new simple argument to derive the vanishing of the 0(1) part of the prefactor. 
According to it the O(l) part of the prefactor vanishes because it is proportional to the determinant of 
the Hessian of the fermionic variables which must be zero due to SUSY breaking. This argument is very 
intuitive and illustrates the similarities between the breaking of a standard bosonic symmetry, that produces 
Goldstone modes and infinite fluctuations, and the breaking of a SUSY, that produces a zero prefactor. 
The above theoretical predictions on the vanishing of the prefactor on the BM solution are also verified 
numerically at the end of section [3] 

Recently Aspelmaier, Bray and Moore have shown that the spectrum of the Hessian of the TAP 
free energy contains an isolated eigenvalue besides the continuous band. They found numerically that this 
eigenvalue vanishes on the BM solution. The presence of a zero eigenvalue in the BM solution is very 
important since it is a possible explanation for the vanishing of the prefactor and for the apparent violation 
of the Morse theorem |13l IH1 1??], see discussion in the conclusions. In sectional we prove rigorously the 
vanishing of the isolated eigenvalue showing that it is another consequence of SUSY breaking. Furthermore 
we discuss the behavior of the eigenvector corresponding to the isolated eigenvalue at the lower band edge 
of the complexity where the continuous band of eigenvalues extends down to zero. At the end we briefly 
discuss the results. 

2 Complexity and Supersymmetry 

We are interested in the solutions of the TAP equations diFTAp{m) = 0, where FTAp{m) is a given model- 
dependent TAP free energy. In particular we want to compute the density of TAP solutions of a given free 
energy, 



a=l 

In the previous expression TAP solutions are labeled by the index a, and {m a } indicates the corresponding 
set of local magnetizations. The density can be expressed as an integral over the whole m-space of a delta 




(1) 
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function of the TAP equations: 



P(f) = E / I[dm z 6(mi - mf) 8[F TAP {m a ) - Nf] = 

]\dm t 6(d z F TAP (m)) \det{d t d 3 F TAP {m))\ S[F TAP (m) - N /] 



(2) 



In order to use an exponential representation of the determinant its modulus is dropped. This is a crucial 
step that corresponds to count each solution with a weight proportional to the sign of its determinant. The 
resulting object is not the number of solutions with free energy / but rather the number of solutions with 
positive determinant minus the number of solutions with negative determinant. Therefore, as we said in 
the introduction, the problem splits in two: from one hand we compute the sum without the modulus but 
on the other hand we want to obtain information on the actual the complexity, i.e. on the sum with the 
modulus. In particular we want to compute the total complexity but we want also to recover the trivial 
result that the sum without the modulus of all solutions be equal to 1 because of the Morse theorem i.e. 
J dfp(f) = 1. We also use an exponential representation for the delta function, 



TAP) — 



i 

det(didjF TAP ) 



J ioo IT^ eiq> nTxidiFTApim)^ , 



dipi dipi exp 



^ tpitpjd t djF TAP (m) 



(3) 
(4) 



where {ip,ip} are anti-commuting Grassmann variables. The delta function over the free energy is also 
expressed as an exponential 



Nf) 



r +i °° du 

/ fi exp[u(F TA p-iV/)] , 



5{F TA p 
In this way we can write .12., 

p{f)= J duVmVxV^V^ eS(m,*,$,i>,«) ; 

where the action S is given by, 

S(m, x, V>, u) = ^XidiF TA p(m) + ^'ip i if} : jdid j F T Ap(m) + u[F TA p{rn)] - uf 



(5) 



(6) 



(7) 



A key property of action Q is its invariance under a generalization ^ [SI E] of the Becchi-Rouet-Stora- 
Tyutin (BRST) supersymmetry [111 I12| : if e is an infinitesimal Grassmann parameter, it is easy to verify 
that the action (JJJ is invariant under the following transformation, 



Srrii = eipi Sxi — —euipi Sipi = —e Xi 
The TAP free energy of the SK model is given by |2] , 

F TA p(m) = -^^Jijmimj + ~y"]<fto(g,mj) , 

ij i 

with, 



(/)o(q,m) 



-(1 +m) log 



-(1 + m) 



-(1 - m) log 



-(1 - to) 



- log(l - to 2 ) + to tanh _1 (m) - log 2 - ^-(1 - q) 2 



(8) 



(9) 



(10) 
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The variables mt are the local magnetizations, and q is the self-overlap of the TAP states, 



q N 

i 

while the quenched couplings are Gaussian random variables with zero mean and variance 1/N. The 
TAP equations and the Hessian of the free energy are respectively, 

f3d l F TAP (m) = -P^Jijrrij + fafami) =0 , (12) 

N 

with 



(3didjF TAP (m) = -0Jij + fo(q,mi) kj — m t m j, (13) 



<j>x{q,m) = 2 (1 - q)m + tanh 1 (m) , (14) 

<h(q,m) = I3 2 {l-q) + —^ . (15) 

1 — m 

The last term of order 0(1/N) in (|13|l is very important for the present discussion. We perform an annealed 
calculation of the number of TAP states, i.e. we directly average p(f) in over the distribution of the 
quenched couplings Jij . Thus, the average number of TAP states becomes, 

pC/J = f VJ P ( J ) Vm Vx V ^ V ^ du e^ m ' x '^' u ^ = J VmVx V4> X></> du e 0S„(m,x,&iM _ ^ 
where 

(3S av = In / -pj p(J) e PS(m,x,M,v.) ( 17 ) 



3 The Macroscopic Action 

The term in the action Q that depends on is 

/3<Sj = — f3 J y - (xiTfij + Xjirii + ipjipi + i>iipj + umiuij) (18) 

By averaging the exponential of Sj with respect to the J's the contribution of the previous term to the 
averaged action S av turns out to be: 

1 — — 

i<j 

— — y^(ximj + Xjirii + ipjipi + tp/ipj + urriimj) 2 - — - ^ (2x i m l + 2-ip^i + um,m,) 2 (19) 
4iv ' — ' J 4iv * — ' 

ij i 

Introducing the following shorthand notation for various macroscopic quantities, 



(20) 



we can express the first term in eq. (|19(l as 

N(3 2 (^+uqR+^^x 2 + ^+uee + 6v + V6-^) ( 21 ) 
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This term is order N while the second term in eq. I|19fl is order 1 since it is the sum divided by N of N 
terms that are local in i: 



(3 2 — 1 — 

— ^2(2x i m l + 2ipiipi + umifrii) 2 = — Fi oc (mj, Xj, j)j , tpj) (22) 



Looking at JJJ, we see that the averaged action l|17|l is given by 

i i i 

+ ^^2Fi oc (m i ,x i ,lp i ,ipi)-20 2 NSe (23) 

i 

The last term comes from the term of order 0(1/ N) in the Hessian i|13|) . Indeed from Q we see that the 
contribution to S av of this term is 

2B 2 — - 

-^rJl ^i^ m 3 m * = -2N(3 2 ee (24) 

ij 

Note that the action H23|) is still invariant under the BRST transformation 

5m l = eifii ; Sx, = -euip, ; Stpi = -ex t ; (25) 

To proceed further we eliminate the squares of the macroscopic quantities R,T,q, 9, 9, V, is, in the action (|23|l . 
In particular the exponentials exp[i? 2 ] and exp[— T 2 ] are expressed as gaussian integrals over respectively 
the real and the complex axes, introducing the variable r and t. We also use an integral representation 
of a delta function on q that introduces an extra parameter A conjugated with q. The fermionic part 
is simplified through the following odd Hubbard-Stratonovich transformation that introduces four new 
parameters {/?, p, ~p, p} 



cxp[— ( 9 f)L^^^] = detL J dpdpdpdp exp[( p p ) L 1 



(26) 



Where according to eq. ll'.'il the matrix L is given by 



L =-" 2 ("i 2 i) ^ 

By means of the various transformations the total action can be expressed as an integral over the variable 
u and eight macroscopic bosonic and fermionic variables = {r,t,<7,A,p,p,/I,/i} 



Where 



J du dQ expLYXl! + 7VE 2 ] (28) 



Si = -/3uf - Xq + qru ^ ,, 



u(3 2 qu(3 2 q 2 uf3 2 q 2 u 2 (3 2 
2f3 2 ' 2/3 2 ~ + ~ 4~ + 4 



+ ]Ifi + (2 - u)pp + pp (29) 
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and 



log 



dmdxdipdip exp x<p(q, m) + ipipfoiq, m) + u<f>o(q, m)- 



qp3 2 x 2 1 — — 2 u fl 2 \2 

— 2 h Fi oc {m, x, ip, ip) + rmx + tipip + Am +—^-{l — q) + 

pplmip — ipm~p(3 — p(3xip — ipx~p/3 



(30) 



Following Kurchan |KS| we note that the original invariance of the action under the BRST transformation 
reflects itself in the fact the macroscopic action Si + £2 is invariant under the following SUSY transforma- 
tion: 



6fi = 


( ru 


up) + qup3+ 2 ± 


5p = 




- | - QU^j e 


8~p = 







5-p = 







Sq = 




2JL 

13 


S\ = 







5r = 




(2pp-pp)e 


St = 


((« 


-2)p3U-pJp)e 



(31) 



To obtain the previous transformation we can consider a generical SUSY transformation and determine 
its parameter in order that £2 be invariant. To do this we must use the fact that £2 is invariant if we 
transform the integrand in 1)30(1 with a transformation like JHj. This happens because of the vanishing of 
the surface terms at m = ±1 and x = ±«oo, as can be easily checked. Thus the variation of £2 induced by 
the transformation (jSJ) on its integrand is a quantity, depending on the macroscopic variables, that must 
sum up to zero and can be added to the variation induced by a generical SUSY transformation between 
the macroscopic variable in order to determine (|31l) . 



4 The Prefactor of the Exponential (I) 

In this section we discuss the general problem of evaluating a SUSY action by means of the saddle point 
method. We are interested in object of this kind: 

J d{$}d{*}exp[7VF($,^)] (32) 

Where {<&} and {"J} are two sets of respectively bosonic and fermionic variables. Expanding F around a 
saddle point with the substitution $ = $0 + 5$/VN and * = * + Sty/Vis we obtain 



NF{*, 9) = NF(* 0l * ) + ^f*f* + 3=^«*f* + O(^) , 03) 
and the integral 1(32(1 over the action can be rewritten as 

exp[7VF($ , *„] J d{$}dm exp[^j^8$8$ + 3=^**** + 0("^)] (34) 

The second factor in the above expression gives the corrections to the leading exponential contribution 
of the saddle point. We note that because of their very nature the integral on the fermionic variables in 
((32(1 in any case does not give a term of the form e^; this feature can be accounted for by saying that the 
fcrmions must be "set to zero" on the saddle point, that is we always have ^0 = 0. 

In the following section we will report a general argument |13| in order to prove that if the action F 
possess a SUSY the prefactor of the exponential contribution on a non-SUSY saddle point is zero at all 
order in an expansion in powers of l/N. In this section instead we provide a simple argument to prove the 
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vanishing of the first term in the expansion. According to the previous expression this term is simply given 
by: 



f n2 p fp, p / a2p \ / f) 2 F \ 

/ d{$}d\m expf^f— <5$(5$ + -= — = det -= — det -f— 



(35) 



The invariancc under the BRST transformation reads 



dF 

* = = DF = for «5$ 4 = T lk & k e (36) 

where «& represents the total set of fermionic and bosonic fields. Deriving the previous expression with 
respect to a generic component of 3> we have 

-??—T tl *, + §-T tm = (37) 

On a saddle point the second term in the l.h.s. of the previous equation is zero by definition therefore, 
if Tm&i is different from zero, it is an eigenvector of the matrix dF/d&kd& m with zero eigenvalue. The 
vector Tki<&i is strictly zero by definition on a SUSY saddle point while it is non-zero on a non-SUSY 
saddle point. The fermionic component of <&, i.e. {71, p} are "set to zero" on both the SUSY and the 
non-SUSY saddle point, as a consequence the bosonic component of the variation vector Tki<&i is strictly 
zero on both saddle point. Thus only the fermionic component of Tki&i on the non-SUSY saddle point 
can be non-zero, this means that the part of the Hessian with a zero eigenvalue is the fermionic one, as a 
consequence: 

det ( -L-) = (38) 

\d*d*J NO-SUSY 

Coming back to equation (|35|l we see that the O(l) contribution to the prefactor of the exponential term 
vanishes on a non-SUSY saddle point. The same argument does not apply to a SUSY saddle point where 
the Tu&i is strictly zero and equation (j37(l is verified trivially. 

In order to verify directly the previous argument we have computed the fermionic determinant expanding 
F = Si + £2 in powers of ^} = {7J, /ip, p}. The resulting expression is 



d 2 F _ ( 1 + f3 2 (i(j^mx) (3 2 (^7 



dmd^ \ 2 - u + /3 2 (tpipx 2 ) 1 + ft 2 (ipipmx) 



(39) 



Where the square brackets mean averages computed with respect to the integrand of £2 feauation l30fl with 
the fermionic variables {jl, p>~p, p} set to zero. Performing the integration over ip,ip and over x we obtain 
the components of the Fermionic Hessian as averages with respect to the variable m. Skipping from the 
macroscopic variable {r, t, A, q, u} to the BM variables {B, A, A, q, u} of [7] we obtain: 



d 2 F f 1 + /3 2 ((1 - m 2 )m mA ~^ h m ) (3 2 ((1 - m 2 )m 2 ) 

Wm> = \ 2-u + f3 2 ((l-m 2 ) {t!iah ~ 1 " l -™ i A)2 - q02 ) l + /3 2 ((l-m 2 )m mA -^ h 



(40) 



Where the square brakets represent averages computed with respect to the standard BM action. Wc checked 
numerically that when computed on the BM solution the determinant of the previous matrix is zero. 



5 The Prefactor of the Exponential (II) 

In this section we briefly recall the argument of Kurchan |13| in order to show that the prefactor vanishes 
at all orders in the expansion in powers of 1/N. In the global variables {#} of the previous section the 
integral over the action reads: 

e JVF(* ) j exp[F 2 ** + — Lf 3 *** + — F 4 ***# + -..], (41) 

J ' V N N 
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and the expansion in powers of l/y/N is: 



e NF(* 0) I ( ) Ar ($) . (42) 

r=0 



Now we make the change of variables 3? = 3?o + «&'/%/ N. The integrand must be invariant under the BRST 
transformation 5&i = eTik&k, that in the new variables reads 



D - = T*#< 0) — + — T,**t — 



^**i 0) 3§7 + T **k-^r = + D * ( 43 ) 



The invariance reads 



— A r (&) = VND A + Y, -7= Po4+i+M) = (44) 



r— v v 7 r— v v 

Since it must be satisfied for each value of N each term in the r.h.s must vanish separately leading to: 

D A = 0; D Q A r+1 + D x A r = (45) 

The operators D and D\ satisfy 

D Q D X + L>iL>o = D 2 Q = D\ = Q (46) 

Note that the equation DqAq = is trivially satisfied on a BRST-susy saddle point, indeed according to eq. 
(I43|l Do depends on the vector of the BRST transformation Tifc$° that is zero by definition on a symmetric 
saddle point leading to D = 0. Instead if D a =/= we may write DqB = 1 with B a constant times one 
of the fermionic variables for which T^^i is non-zero. Multiplying by (D B) k — 1 and integrating 
repeatedly by parts we have 



J d{*'}4.(*') - J d{&}(DoB) r+1 A r (<S>') = 
= J di&jiDoByA^^B) 



d{&}D BA {DxB) r = 
J di&jBDMD^y = (47) 



Thus looking to l|42|l and H47JI we see that on a non-SUSY saddle point the prefactor of the exponential 
vanishes at all order in powers of l/N. As noted in J3| this may happens either if the prefactor is strictly 
zero at all orders or if it is of the form e~ aN with a > 0. 

Let us note that the action (|28|l is not precisely of the form l|32() since £2 defined in (|30|l contains an 
extra dependence on l/N coming from the local terms (|22|l . These terms obviously enter the l/N expansion 
but they do not change the result on the vanishing of the expansion at all orders. Indeed the presence of 
these extra terms in (I41f) does not change the form (|42|) and the subsequent derivation. The same is true 
for the argument of the previous section, indeed these terms give a multiplicative correction different from 
one to the prefactor that cannot change the final result that depends on the vanishing of the fermionic 
determinant. However it is important to note that these terms should be considered in the construction of 
l/N corrections to the real complexity, i.e. the one with the modulus. 



6 The Isolated Eigenvalue 

The density of eigenvalues corresponding to the Hessian without the last term in (|13fl can be obtained 
through standard theorems on random matrices (see |151 \lb\). It turns out to have always a positive 
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support with a lower band edge going to zero if the quantity 



= 1-/3 2 £(1-^) 2 



goes to zero. The BM solution has a finite value of xp meaning that the corresponding TAP solutions has 
a definite positive Hessian. On the other hand the Morse theorem predicts that there is an equal number 
of solutions with positive and negative determinant of the Hessian, so the problem is: where are these 
negative determinant solutions? Note that the answer cannot come from the finite- TV tail of the eigenvalue 
density if we assume that at finite N the eigenvalues outside the support are concentrated at a distance 
f rom lower band edge, as in most cases in random matrix theory |18|. 

In JU] Aspelmaier, Bray and Moore (ABM) has shown that the inclusion of the last term in the Hessian 
(|13fl modifies the continuous spectrum of the first two terms with the addition of an isolated eigenvalue. 
Most importantly they checked numerically that this isolated eigenvalue vanishes on the BM solution solving 
the apparent inconsistency with the Morse theorem. In this section we prove rigorously the vanishing of 
the isolated eigenvalue showing that it is a consequence of SUSY breaking. 

The appearance of the isolated eigenvalue is connected to the fact that the last term in the Hessian 1)13(1 
is proportional to a projector. Given two symmetric matrices A and B that differ by a projector P = \a) (a\ 
we have: 

A = B + P — >detA= (1 + (a|B" 1 |a))dctB (48) 

This relation can be easily proved considering an orthonormal base with, say, the first element equal to 
\a), such that the only non-zero component of Pjj is P±i = 1. The same argument applies to B = A — P 
leading to: 

1 + (alB- 1 \a) = - 1 (49) 
1 — (a\A L \a) 

Thus the projector modifies the spectrum, in particular if (1 + (a|i? _1 |a)) is zero the isolated eigenvalue 
vanishes and the corresponding eigenvector is £? _1 |a:). 

The quantity (1 + (a|_B _1 |a)) is essentially the ABM quantity 1 — 20 2 H that must be zero in order 
for the isolated eigenvalue to vanish JHj- By direct inspection one can check that the expression for the 
determinant of the fermionic Hessian eq. (|4L)fl turns out to be proportional to 1 — 2(3 2 H, see equation (8-11) 
in |19) . Since in section 0] we proved that the fermionic determinant must vanish on a non-SUSY saddle 
point like the BM one this also prove rigorously that the isolated eigenvalue vanishes: 

-)2 



/ a2p \ 

1 - 2(3 2 H L_ HTTQV oc det ^ — =0. (50) 

I NO-SUSY U*Wmo_«tt S v 



NO-SUSY 

In the following we prove directly the connection between the isolated eigenvalue and the determinant of 
the fermions. The relation H49[) allows to consider the quantity 1 — (alA -1 |a) instead of (1 + (a\B^ \a)). 
That is we can use the full Hessian 

A i:j = d 2 F/dm l dm :j , (51) 

instead of the Hessian without the projector term. Therefore we must evaluate (m\A~ 1 \m) /N on all 
solutions of the TAP equations: 

12 a -Nl2ij m i(-^k).. m J _ J jy m^m,- exp[5[m, x, j>, g] _ J Mgexpjgg 
J2 a fexp[S[m,x,tp,i)]\ /exp[/3S] 

Where the action f3S is given by eq. i|23|) . The presence of the term 90 in the integrand must be taken 
into account when performing the decoupling of the variables {6,9,u,V}. In general this can be achieved 
by multiplying the integrand in the l.h.s. of eq. for some proper prefactor depending on the matrix L 
and on the integration variables {n,JI, p,~p}. In our case L is given by eq. (|27|l and the prefactor turns out 
to be fip. Thus we obtain that: 

■Hr.,*!^" (53) 

N J c(B exp[iV hi + TVl^aJ 



9 



By computing the previous quantity at the saddle point we obtain that the bosonic contributions of order 
e* 5 " and of order 1 (given by the bosonic determinant, see eq. Q35[l) are equal in the numerator and in 
the denominator. Instead the fermionic contributions are different. In the denominator the contribution is 
given by the fermionic determinant as in (|35|) while in the numerator the presence in the integrand of the 
factor nJI gives only a contribution equal to d 2 F/dp&p that is a finite quantity proportional to the diagonal 
terms of eq. (|4U|) . Thus we have that: 

1 r) 2 F / r) 2 F \ -1 / r) 2 F \ 

— (mlA-V)* det f -== — — ►l + (aLB- 1 |a}ocdet[-^= — (54) 

Going back to equation (|48J) we see that the determinant of the Hessian is proportional to the determinant 
of the fermionic components of the macroscopic action and vanishes on the non-SUSY solutions. 

According to the prediction of the BM solution the determinant of the Hessian TAP solutions have a 
zero eigenvalue at any value of the free energy. When the free energy goes to the equilibrium value the 
continuous band develops itself a zero eigenvalue and it is interesting to ask if the two eigenvectors, the 
one of the continuous part and the one of the isolated eigenvalue, are different. 

Let us reconsider the expression for the determinant of the matrix A defined above in equation (|48|l . 
We rewrite it in a way that takes into account the possibility that det B = 0. 

det A = det B + (m\C\m) , (55) 

where C is the matrix of the cofactors of B such that B^ 1 = Cj dctB. If dctB is different from zero the 
condition for the vanishing of det A is that the l.h.s. of equation Q49|l vanishes and the eigenvalue of A with 
zero eigenvector is B~ x \m). Instead if det B = the condition for the vanishing of det A is simply 

(m\e ) = , (56) 

where \eo) is the eigenvector of B with zero eigenvalue; the previous result can be simply derived considering 
the form of the matrix C in a basis in which B is diagonal. The last equation implies that the eigenvector 
|eo) is also an eigenvector of A with zero eigenvalue. Instead the vector C\m) is simply zero. To understand 
what is the behavior of the eigenvector of the isolated eigenvalue it is convenient to study the situation in 
which the matrix B has a small eigenvalue equal to e corresponding to an eigenvector |e e ). In this case one 
finds that the eigenvector of the matrix A corresponding to the vanishing isolated eigenvector is equal to 

B _1 |m> oc|e £ )+0(Vi) . (57) 

As an eigenvalue of B approaches zero the eigenvector i3 -1 |m) of A becomes equal to the eigenvector of 
B with zero eigenvalue. The previous arguments can be extended in the case when there are more than 
one zero eigenvalue in the matrix B. If the matrix B has k eigenvectors | e^} with zero eigenvalues the 
matrix A has at least k — 1 eigenvectors with zero eigenvalue. Instead if the condition (|56|l holds for all 
the k eigenvectors the matrix A has k eigenvectors with zero eigenvalue. Finally let us note that the 
condition l|56(l implies that the eigenvector with zero eigenvalues of the TAP solutions corresponding to the 
equilibrium states is orthonormal to the vector of the magnetizations. 



7 Conclusions 

In the introduction we saw that by dropping the modulus of the determinant we face two problems. 
Assuming that the BM theory gives the actual complexity (the one with the modulus) the problem is to 
obtain the exact result for the sum of solutions weighted with the sign of the determinant; for instance 
the sum over all solutions must yield 1 because of the Morse theorem. In ^3] Kurchan suggested that 
this is recovered in the BM theory due to the vanishing of the prefactor at all order in 1/N. Indeed if the 
prefactor is proportional to e~ NT,BM (such that e ~~ NT,BM e NT,BM — 1) the resulting expansion in powers of 
1/N is zero at all orders. According to our results the prefactor vanishes at all orders not only for the 
total number of solutions but separately at any given free energy /, therefore, applying the same argument, 
we must conclude that at any given free energy there is an equal number of solutions with positive and 
negative determinant, much as if the Morse theorem be valid at any / separately. 
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The previous arguments imply that the BM saddle point describes solutions with both positive and 
negative Hessian at all free energies. On the other hand the leading part of the spectrum of the TAP Hessian 
is strictly positive in the BM theory and the question is: where are the negative Hessian solutions? Recently 
Aspelmeier, Bray and Moore ^5] have shown that a solution to this problem can be found noticing that 
the TAP spectrum contains a single isolated eigenvalue besides the continuous strictly positive band. They 
have checked numerically that this eigenvalue is zero in the BM theory. Therefore in the thermodynamic 
limit the Morse theorem is recovered because the solutions have zero determinant. At finite N the isolated 
eigenvalue acquires a non-zero value 0(l/\/~N), each solutions becomes a minimum or a saddle and at 
the same time another solution with an opposite value of the isolated eigenvalue develops at a distance 
0(1/ \/N). These results are valid for the total number of solutions (consistently with the Morse theorem) 
but also with the solutions at any given free energy consistently with the conclusions of the previous 
paragraph. Here we have proved that the vanishing of the isolated eigenvalue at any given free energy is 
an exact result that is a consequence of SUSY breaking. 

Since it has been shown JHj that the SUSY solution is inconsistent (except at the lower-band edge, 
see jHlEIj) in full replica-symmetry- breaking (FRSB) models like the SK model (at variance with 1RSB 
models), the BM solution remains the only viable candidate to describe a finite complexity in these models 
and seems to be supported by recent numerical results [20] ■ The SUSY approach has given us further insight 
into this non-SUSY theory, due to the presence of the isolated eigenvalue founded in J5j we know that the 
TAP solutions described by the BM saddle point are inflection points in the TAP landscape while they split 
in couples minima-saddles at finite N. These features may be of some importance for the dynamics. Note 
that the presence of zero eigenvalues in the spectrum is a well-known feature of the lowest TAP solutions 
but has a different origin, indeed it is due to the fact that at the lower band edge the continuous band of 
eigenvalues extends down to zero. We have considered the problem of the relation between the eigenvector 
associated to the vanishing isolated eigenvalue and the eigenvector associated to the zero eigenvalue coming 
from the continuous band . In the end of section JHl we have shown that the two eigenvectors concide. 
However in this case the SUSY approach yields the non-trivial result that the eigenvector of the Hessian 
corresponding to the zero eigenvalue is orthonormal to the vector ot the magnetization {m}, according to 
equation 

An open problem of the BM theory is how to compute the 1/N correction to the number of solutions. 
A possible way is to remove the effect of the vanishing isolated eigenvaue responsible for the cancellation 
of the prefactor. Work is in progress in this direction and on the extension of the results presented here to 
the replicated case. 
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